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ABSTRACT 

Using photometric data from the 2MASS second incremental release point 
source catalog we constructed Kg vs. (J—Ks) color-magnitude diagrams (CMDs) 
of Galactic globular clusters (GGCs) for which the JKg photometric data have 
been made available up to now. On the CMDs of 13 GGCs we identified RGB 
bump features and derived luminosities of bumps in Ks {Kf^'^^) and in Mk, 
(M^"™'') for 11 of them. We reconfirm the relation between M^"™'' and metal- 
licity [Fe/H]cG97 or [M/H] such that the luminosity of bump becomes brighter as 
metallicity [Fe/H]cG97 or [M/H] decreases. This result is similar to that obtained 
by Ferraro et al. (1999) for the relation between M^"™^ and [M/H] based on 
observations of 47 GGCs which were conducted in the optical region (in V vs. 
(B — V) CMDs). Our results show the same trend as those found by Ferraro et 
al. (2000) for the relations between M^*^™^ and metallicity [Fe/H]cg97 and [M/H] 
derived from observations of 8 GGCs in the same near-infrared region {K vs. 
(J — K) CMDs). Combining with the data of Ferraro et al. (2000), we derive a 
robust relation for the metal- dependent luminosities of the bumps. 

Subject headings: globular clusters: general — stars: color-magnitude diagrams — 
stars: luminosity function — stars: evolution — stars: population II 



1. Introduction 

The red giant branch (RGB) bump in globular clusters was theoretically predicted by 
Thomas (1967) and Iben (1968) as a region where evolution through the RGB is stalled for 
a time when the H burning shell passes the H abundance inhomogeneity envelope. This is 
produced by the stellar outer convection zone at the H shell burning RGB stage after the 
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first dredge-up of a star in a globular cluster. The first convincing identification of the bump 
was that of the metal-rich cluster 47 Tuc (King, Da Costa, & Dcmarque 1985). The position 
in luminosity of the RGB bump is a function of metal abundance, helium abundance, and 
stellar mass (and hence cluster age) as well as any additional parameters that determine the 
maximum inward extent of the convection envelope or the position of the H burning shell. 
Thus the position of this peak on the giant branch of GGCs should provide observational 
constraints on these parameters. 

RGB bumps in Galactic globular clusters (GGCs) were first obscrvationally identified 
in a systematic way by Fusi Pccci ct al. (1990). They found RGB bumps in 11 GGCs 
from the peaks in the differential luminosity function and a change in the slope of the 
integrated luminosity function using V vs. {B — V) color-magnitude diagrams (CMDs). 
They discovered a correlation between the positions of RGB bumps in the GGCs and their 
metallicities. Subsequent to this work, many studies on this subject have been carried out 
(Sarajedini & Norris 1994; Brocato et al. 1996; Savianc ct al. 1998). Recently, Ferraro et al. 
(1999) systematically identified the positions of RGB bumps in 47 GGCs in V vs. {B — V) 
CMDs and confirmed the strong correlation between the absolute V magnitudes of RGB 
bumps and metallicities of GGCs. 

In the near-infrared band, the the locations of RGB bumps were identified for 8 GGCs 
in K vs. (J — K) CMDs and a similar correlation between absolute K magnitudes of RGB 
bumps and metallicities was derived by Ferraro et al. (2000). 

The advantage of observing GGCs in the near-IR is an enormous reduction in reddening. 
The interstellar reddening in the K band (Ak) is roughly 10% of that in the visual {Ay) 
part of the spectrum. This imparts a great advantage to near-infrared observations especially 
toward the disk and severely obscured bulge directions. Therefore, near-infrared observations 
are much preferred for disk and severely obscured bulge GGCs and for as yet unknown GGCs 
severely reddened by the Galactic bulge in the visual bands. Near-infrared observations are 
also useful for normal disk and halo GGCs because the {V — K) color index is a good 
temperature indicator. 

In the near-infrared band, RGB bumps of GGCs can be a very useful distance indicator 
because they are relatively bright, and their position in CMDs is close to the HB (Horizontal 
Branch). Moreover, the HB, which is a good distance indicator in the V band is not useful in 
the near-infrared CMDs. In the near-IR CMDs, the HB is diagonally slanted from the upper 
right to the lower left position and is not horizontal at all. So it is difficult to measure the 
HB levels of GGCs in the near-IR CMDs. It is expected that the RGB bump is a distance 
indicator which would have the same amount of observational error as the HB used in the 
optical region since it has a similar HB metallicity dependency. 
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The goal of this work is to obtain the absolute magnitudes of RGB bumps in as many 
GGCs as possible in the near-IR band and to derive a robust relation between the luminosities 
of the bumps and metallicities of GGCs for use as distance indicators in the IR bands. Using 
the photometric data from the 2MASS (Two Micron All Sky Survey) second incremental 
release point source catalog we have investigated the existence of RGB bumps in GGCs and 
have found a correlation between the absolute magnitudes of RGB bumps and metallicities 
of GGCs, comparing our findings with the results of Ferraro et al. (2000). Combining with 
their data we provide a more robust correlation for the use of RGB bumps as standard 
candles for deriving distances to other GGCs. 

In § 2 we discuss some characteristics of GGCs' near-infrared CMDs and in § 3 investi- 
gate the location in luminosity of the RGB bumps and the relation between the luminosity 
of the bump and metallicity of GGCs, and in § 4 we briefly summarize our results. 

2. IR Color-Magnitude Diagrams 

Wc used photometric data from the 2MASS second incremental release point source cat- 
alog which were obtained at two 1.3-m telescopes in the northern and southern hemispheres 
in the near-infrared bands at J (1.25/im), H (1.65/im), and Ks (2.17/xm; K short = medium 
modifled K). Limiting magnitudes of the 2MASS second incremental release point source 
catalog photometric data are about 15 mag and the released data from 2MASS now covers 
~47% of the sky. 

The JKg photometric data for GGCs in Harris' (1996) catalog which are relatively close 
to the Sun are obtained from the IPAC (Infrared Processing and Analysis Center) and Kg 
vs. {J — Kg) CMDs are constructed. We rejected data whose errors are larger than 0.15 mag. 
We can construct CMDs of about 46.8% of GGCs in Harris' (1996) catalog from the 2MASS 
second incremental release point source catalog. However, the quality of many CMDs is not 
good enough to find the bumps because the limiting magnitudes are too bright. 

In Figure 1 we present CMDs of 11 GGCs for which the RGB bump positions are 
accurately measured by the systematic analytical method described in § 3.1. According to 
Figure 1 the major characteristics of the CMDs of GGCs in the Kg vs. (J — Kg) plane are 
as follows. 

First, the BHB (Blue Horizontal Branch) and RHB (Red Horizontal Branch) arc diag- 
onally tilted from the upper right to the lower left position and such features are obvious 
in the CMDs of M4 (NGC 6121) and M107 (NGC 6171) which have both BHB and RHB. 
Even the RHB of 47 Tuc (NGC 104), which has only RHB, is clearly tilted. So in the Kg vs. 
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{J — Ks) CMDs it is difficult to determine the HB level of any GGC and to correlate RGB 
bump positions and HB levels in the near-infrared CMDs, in contrast to the case for optical 
CMDs. 

Second, the brightness interval between the MSTO (Main Sequence Turnoff) and RGB 
tip is larger (~9 mag) than in optical CMDs (~6.5 mag), as can be easily seen in M4. So the 
magnitude resolution in the Kg vs. {J — Kg) CMDs is larger than in optical CMDs. However, 
color resolution in the Kg vs. {J — Kg) CMDs is smaller than in optical CMDs because in 
the latter the RGB and AGB (Asymptotic Giant Branch) are separated from one another 
at least in the lower part of the AGB. In contrast, in the Kg vs. (J — Kg) CMDs the RGB 
and AGB overlap, as clearly shown in the CMD of 47 Tuc whose AGB is relatively rich. In 
the optical CMDs of NGC 362, 47 Tuc, and M71 (NGC 6838), the separation of their RGBs 
and RHBs is clearly seen but they partially overlap in the Kg vs. (J — Kg) CMDs. This also 
results from the fact that the color resolution in the Kg vs. (J — Kg) CMDs is lower than 
that in optical CMDs. 

Third, in the case of M22 (NGC 6656), its RGB is broader than those of other GGCs 
as found in the optical CMD (Peterson & Cudworth 1994). Moreover, in the lower part of 
the RGB, contamination from bulge component stars is severe. To the right of the RGB of 

M22 there exists another RGB component. This RGB component would be due to bulge 
RGB stars, which extend prominently rcdward in the optical CMD (in our unpublished data) 
because of the strong blanketing effects of heavy metals as found in metal-rich GGCs such 
as NGC 6553 (Ortolani, Barbuy, & Bica 1990). 

Last, because NGC 362 and 47 Tuc CMDs are located in the direction of the SMC, a 
contribution from the SMC appears in the lower right region of their CMDs. 

3. The RGB Bump 

3.1. The Luminosity of the RGB Bump 

In order to accurately measure the luminosity of the RGB bump and construct lumi- 
nosity functions of the GGCs, we applied several standard procedures to delineate only the 
RGB sequences for all GGCs except for M22 and M71. First we rejected visually clear HB 
stars and AGB stars in the process of distinguishing outlying field stars from RGB stars. 
Second, by binning the RGB sequences in 0.5 mag intervals we measured the average and 
sigma of each bin and rejected stars 2(t away from the mean value of each bin. We then 
remeasured the average and sigma of each bin and further rejected stars 2(7 away from the 
remeasured average of each bin, and iterated this process until no stars were rejected and 
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the average and sigma of each bin converged. Third, from each RGB sequence after field 
star rejection, we constructed a differential luminosity function and an integral luminosity 
function for each GGC following the classical method of Fusi Pecci et al. (1990). 

Then, we searched for a significant peak in the differential luminosity function and a 
corresponding large slope change in the integral luminosity function simultaneously for each 
GGC. Where the significant peak position in the differential luminosity function and the 
position of large slope change in the integral luminosity function coincide, we identified the 
significant peak as the RGB bump of the given GGC and measured this position {Kf™'^^). 
However, in several GGCs the existence of RGB bumps is very clear only when considering 
the differential luminosity functions. 

In the case of M22, since its CMD is severely contaminated by field and bulge stars, we 
clipped off the proper region in order to isolate the RGB sequence to Kg = 13. Below Kg 
= 13, isolation of the RGB sequence is difficult because of severe contamination by the field 
and bulge stars. 

In the case of M71, the RGB is too poorly populated to be clearly defined. So in order 
to isolate the RGB sequence of M71 we used the 2a clipped RGB sequence of 47 Tuc with 
a similar metallicity as a template for the RGB sequence of M71. The remaining processes 
were the same as for the other GGCs. 

Differential luminosity functions and integral luminosity functions for each GGC are 
shown in Figure 2. 

In order to derive absolute magnitudes of GGC RGB bumps (A/]^^™^) from their appar- 
ent magnitudes we need absolute distance moduli and interstellar reddenings. We take these 
values from Table 2 of Fcrraro et al. (1999) except for M22 and M2 (NGC 7089). Ferraro 
et al. (1999) derived absolute distance moduli in a very systematic way using the HB level 
as a distance indicator, or more strictly speaking the ZAHB (Zero Age Horizontal Branch) 
level. They derived ZAHB levels for GGCs by comparing the CMDs of GGCs with synthetic 
H-R diagrams. Then they derived absolute distance moduli using absolute magnitude levels 
of the ZAHB (My^^^) and interstellar reddening. For interstellar reddening they adopted 
values from the compilation of Harris (1996). 

In the cases of M22 and M2, which are not in the hst of Ferraro et al. (1999), we 
calculated their absolute distance moduli based on the same method described in Ferraro et 
al. (1999), referring to other sources for the physical values. For derivation of the absolute 
distance moduli of M22 and M2 we used the magnitude difference between the apparent 
visual magnitude of the ZAHB (Vzahb) and the absolute visual magnitude of the ZAHB 
^^ZAHB^ of each cluster. That is as follows: 
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(m - M)C^^97 ^ _ -Ay^ (FzAHB - M^^™) - Ay. (1) 

For ^zAHB we used equation (2) in § 5 of Ferraro et al. (1999): 

^ZAHB = <Vhb> + 0.106[M/H]^ + 0.236[M/H] + 0.193. (2) 

By the way, a-element enhanced global metaUicity [M/H] values for M22 and M2 were 
computed by equation (1) in § 3.4 of Ferraro et al. (1999): 

[M/H] = [Fe/H] + log(0.638/a + 0.362), (3) 

where fa is the enhancement factor of the a-elements. According to Ferraro et al. 

(1999), for the 19 GGCs with [a/Fe] hsted by Salaris k Cassisi (1996) or Carney (1996), 
they adopted /„ = lO'^/^'^l. For all the others they assumed /„ = lO^-^^ if [Fe/H] < -0.8 
and fa = IQ-o-nFo/u] if [Fe/R] > -0.8. So, for M22 we adopted [a/Fe] = 0.32 from Carney 
(1996), averaged over 9 stars which include all 3 stars of M22 in Salaris & Cassisi (1996). 
For M2 [a/Fe] is not listed in Salaris & Cassisi (1996) or Carney (1996). So, we adopted fa 
= 10°-28 because [Fe/H] of M2 is -1.46. For M^^^^ we used equation (4) for [Fe/H] in § 6.4 
of Ferraro et al. (1999): 

^ZAHB _ ^_QQQ5 ^ 0.3485[Fe/H] + 0.0458[Fe/H]^ (4) 

For Ay we adopted Ay = 3.08E{B - V) from He et al. (1995) and for Ak,, which is 
used in the calculation of apparent Kg magnitude distance moduli for program GGCs, we 
adopted Ak, = 0.32E{B - V) from He et al. (1995). All parameters of M22 and M2 are 
listed in Table 1. 

All these final values for absolute distance moduh and interstellar reddening of program 
GGCs are listed in Table 2. According to Ferraro et al. (1999) the global uncertainty of 
the absolute distance moduli is of the order of about 0.2 mag considering that the derived 
absolute distance moduli are affected by many uncertainties (namely, the evaluation of the 
ZAHB level, the zero point and dependence on mctallicity of the ZAHB level, reddening, 
etc.). Finally, using absolute distance moduh and interstellar reddenings we derived absolute 
magnitudes of the RGB bumps of GGCs (M^""P) hsting them in Table 2 together with 
^Bump Table 2, errors in column (6) are measurement errors and errors in column (7) 
are a combination of measurement errors in column (6) and the global uncertainty of the 
absolute distance moduli, which is equal to 0.2 mag. 
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In the case of M2, the RGB bump appears in the expected position according to equa- 
tions (5a), (5b) and (6a), (6b) in the next section, but the bump is relatively broader com- 
pared to those of other GGCs, and the quality of its IR CMD is poor compared to the optical 
CMD of Lee & Carney (1999). 

Cudworth (1985) claimed that in the optical CMD of M71 there is a strong clump just 
below the HB whose position in the V vs. {B — V) plane is ~ 15 and (B — V) ^ 1.2. 
Transforming these values into the J vs. (J — Kg) plane by the color transformation table 
of BesseU, CasteUi, & Plez (1998), V ^ 15 transforms into Ks ^11.92 and {B -V) ^ 1.2 
transforms into (J — Kg) ~ 0.77. Also, according to Girardi et al.'s (2000) evolutionary 
tracks the former values correspond to ~ 11.94 and {J — Kg) ~ 0.72. These transformed 
Kg and (J — Kg) values of the strong clump of M71 in the V vs. {B — V) CMD nearly match 
the RGB bump position Kg — 11.95 and {J — Kg) — 0.72 derived in this work. So the strong 
clump reported by Cudworth (1985) in the V vs. {B - V) CMD of M71 must be the RGB 
bump derived in this work in the Kg vs. (J — Kg) CMD. 

In the cases of M69 (NGC 6637) and u; Gen (NGC 5139) we can see the RGB bump 
features in their CMDs. However, for M69 it is difficult to measure the RGB bump position 
by systematic analysis because it exhibits a CMD which is too scatted to analyze. Also for 

(J Gen we abandoned the analysis because uj Gen has a wide RGB suggesting a metallicity 
spread (Lcc ct al. 1999; Pancino et al. 2000), which renders measurement of RGB bump 
position meaningless in this work. 

However, in the cases of M30 (NGC 7099) and M55 (NGC 6809), we cannot detect RGB 
bumps at the expected positions taking into account their metallicity after careful analysis 
of their differential and integral luminosity functions. Figure 3 shows their Kg vs. (J — Kg) 
CMDs, differential luminosity functions, and integral luminosity functions. For M30 there 
are hints of bumps at Kg ~ 12.5 and ~ 13.5. However, no slope changes around those 
magnitudes are found in the integral luminosity function. For M55 neither a bump nor a 
significant slope change are found. 

3.2. The Relation between the luminosity of RGB bump and metallicity 

There are two widely used metallicity scales, the Zinn & West scale (Zinn & West 1984; 
Zinn 1980, 1985; Armandroff & Zinn 1988) and the Carretta & Gratton scale (Carretta & 
Gratton 1997). The Zinn & West scale (hereafter ZW scale) employs the most complete 
data set but is not based on high resolution spectra. However, the Carretta & Gratton scale 
(hereafter CG97 scale) is based on systematic high resolution spectra obtained by their own 
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team, and there is a transformation relation between the ZW scale and the CG97 scale in 
equation (7) of Carretta & Gratton (1997). The CG97 scale is more robust since it relies 
on recent high resolution spectroscopic measurements and up-to-date atmospheric models. 
Moreover, in order to compare our results directly with those of Ferraro et al. (2000) which 
employed the CG97 scale, we adopted the CG97 metallicity scale. 

The CG97 scale metallicity of each GGC is given in Table 2. Since only 24 GGCs have 
metallicity determinations from direct high resolution spectra in the CG97 scale, the metal- 
licities of NGC 1851 and M107 are interpolated from equation (7) of Carretta & Gratton 
(1997), and the metallicity of M2 is adopted from the value determined by the morphological 
parameters of RGB stars by Lee & Carney (1999). Recently, Carretta & Gratton's group 
pubhshed high resolution spectroscopic metaUicities of NGC 6553 (Cohen et al. 1999; Car- 
retta et al 2001) and NGC 6528 (Carretta et al. 2001, hereafter CG), extending the CG97 
metaUicity scale into the high metaUicity region ([Fe/H]cG ~ 0.00), and established a new 
transformation equation between ZW scale metallicity and CG scale metallicity. However, 
in the metallicity range of [Fe/H]cG97 = —2.12 ~ —0.70 the interpolated metallicity values 
are coincident within 0.1 dex in the two different transformation equations and it seems that 
transformation equation (7) of Carretta & Gratton (1997) is tighter than transformation 
equation (3) of Carretta et al. (2001). We did not revise the metaUicities of NGC 1851, 
M107, M2, M55, and M69, which have no metallicity measurements from direct high reso- 
lution spectra in the CG97 metallicity scale. But for metaUicities of NGC 6553 and NGC 
6528, which were extrapolated by Ferraro et al. (1999) since their ZW scale metaUicities are 
outside the validity range of the transformation to the CG97 scale metaUicities by equation 
(7) of Carretta & Gratton (1997), we adopted the new high resolution spectroscopic mea- 
surements of Carretta et al. (2001). MetaUicities of the other GGCs are taken from the 
directly determined values of Carretta & Gratton (1997). 

We also consider the global metallicity [M/H] which incorporates a-element enhance- 
ment into the CG97 metallicity scale [Fe/H]cG97 according to equation (1) of Ferraro et al. 
(1999). In the case of [M/H] we directly adopted the values in Table 1 of Ferraro et al. 
(1999) except for M22, M2, NGC 6553, and NGC 6528. In the cases of M22 and M2 we 
calculated [M/H] as described in § 3.1 and in the cases of NGC 6553 and NGC 6528 we first 
calculated [a/Fe] using equation (2) of Carney (1996) with the a-element abundances listed 
in Carretta et al. (2001) and then calculated [M/H] using equation (3) in § 3.1. Values of 
[a/Fe] of NGC 6553 and NGC 6528 are 0.20 ± 0.07 and 0.21 ± 0.04, respectively. 

When we compare M^""'^ with the CG97 metaUicity scale [Fe/H] 

CG97 or the a-element 

enhanced global metaUicity scale [M/H] given in Table 2, we find clear quadratic correla- 
tions as shown in Figure 4. The relations expressed in equations (5a) and (5b) imply that 
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as metallicity decreases, the luminosity of the RGB bump becomes brighter, which is an 
identical result to that of Ferraro et al. (2000). 

^Bump ^ 26 ± 0.16)[Fc/H]cG97 + (1.43 ± 0.45)[Fc/H]cG97 - (0.35 ±0.11) (5a) 
^Bump ^ (^Q_33 ^ o.l8)[M/H]^ + (1.51 ± 0.44)[M/H] - (0.51 ± 0.11) (5b) 

The photometric data of Ferraro et al. (2000) were homogeneous, making use of Glass' 
standard stars (Ferraro et al. 1994; Montegriffo et al. 1995). The 2MASS photometric data 
are also homogeneous and have undergone Global Photometric Calibration (Nikolaev et al. 
2000). Therefore, if we estimate the photometric zero point differences between them we 
can combine the two sets of photometric data by correcting for the zero point difference. 
Comparison of RGB bump positions of clusters common to both works can be used for 
estimation of the zero point difference. However, since the two works were conducted in 
different photometric systems, we first have to render the photometric data in the same 
sustem, that is to say in the 2MASS system. 

We transformed Ferraro et al.'s (2000) results into the 2MASS system using equation 
(Al) of Carpenter (2001). Equation (Al) of Carpenter (2001) can transform the Bessell & 
Brett (1988) system into the 2MASS system. The Bessell & Brett system is the homogenized 
Glass system, and the photometric data of Ferraro et al. (2000) are in the Glass system since 
the data are standardized by Glass' standard stars. So, we transformed the data of Ferraro 
et al. (2000) in the Glass system into the 2MASS system using equation (Al) of Carpenter 
(2001). 

The 3 GGCs of this work overlap with those of Ferraro et al. (2000), and the RGB 
bump luminosities of the two works are listed in Table 3. We list the luminosities of the 
RGB bumps in the Glass system which are original and those in the 2MASS system after 
transformation from the Glass system for the 8 GGCs of Ferraro et al. (2000) in Table 4. 
According to Table 3, in the cases of M15 (NGC 7078) and 47 Tuc, the RGB bump positions 
in this work are fainter by 0.09 mag than those of Ferraro et al. (2000), and in the case of 
Ml 07, the RGB bump position in this work is brighter by 0.01 mag than that of Ferraro et 
al. (2000). However, since these differences are smaller than the errors of the data, we can 
neglect the zero point difference between the two works. Therefore, we combined the data in 
the two works without applying any zero point correction. The results are shown in Figure 
5. 

The sample of 8 GGCs in Ferraro et al. (2000) covers a wide range in metallicity despite 
its small size. While our sample is a little larger, comprising 11 GGCs, it covers a hmited 
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range in metallicity. Therefore, the combined data set increases the sample of GGCs to 16 
and covers a wider range in metallicity than that of the sample GGCs of our work. We 
have found clear relations between GGC metallicity and luminosity of the RGB bump, as 
shown in Figure 5. By taking error-weighted mean averages for the 3 common GGCs and 
combining the remaining data we have the reduced regression equations (6a) and (6b) and 
we plot the results in Figure 5 as solid lines. 

^Bump ^ ^ 0.06)[Fe/H]cG97 + (0.89 ± 0.12)[Fe/H]cG97 - (0.67 ± 0.10) (6a) 
^Bump ^ (0.07 ± 0.06) [M/H]^ + (0.93 ± 0.11) [M/H] - (0.78 ± 0.10) (6b) 

Dashed lines are from equations given in Figure 13 of Ferraro et al. (2000) transformed 
into the 2MASS system and dot-dashed lines are from equations (5a) and (5b). Although 
dashed lines and dot-dashed lines are not exactly coincident with each other, their forms 
are very similar. The first terms of equations (6a) and (6b) are negligible if the errors in 
the parameters are considered. Therefore solid lines are nearly linear and deviate from the 
dashed fines sfigfitly at the metal rich ends. This is due to metallicity revision of NGC 6553 
and NGC 6528. Their original values in Ferraro et al. (2000), which were extrapolated to 
be [Fe/H]cG97 = -0.44 and -0.38, and [M/H] = -0.36 and -0.31, respectively, have been 
revised to [Fe/H]cG97 = -0.06 and 0.07, and [M/H] = 0.08 and 0.22, respectively. 

The overall trend for the RGB bump positions of GGCs to become brighter with de- 
creasing metallicity has been supported by many theoretical models derived from the first 
such model estabfished by Sweigart (1978). However, there is also a moderate helium abun- 
dance dependency (Sweigart 1978) and a weak age dependency (Ferraro et al 1999; Yi et al. 
2001) of RGB bump positions for a given metalficity. 

According to Yi et al. (2001) and Ferraro et al. (1999), the RGB bump luminosity varies 
with metallicity by ^"^^h] ~ 0.96 in the metalficity range [Fe/H] = —2.3 ~ 0.0 and in the age 

range of 7 Gyr ~ 16 Gyr. Yun & Lee (1979) found that RGB bump bolometric luminosity 
varied with [Fe/H] according to a nearly constant relation in the helium abundance range Y 
= 0.1 ~ 0.3. in the theoretical himinosity function analysis of 46 RGB models of Sweigart 
& Gross (1978). Therefore, we assume that the RGB bump luminosity varies with [Fe/H] 
constantly by ^^J"^^ ~ 0.96 in the metallicity range [Fe/H] = —2.3 ~ 0.0, in the age range 
of 7 Gyr ~ 16 Gyr, and in the helium abundance range F = 0.1 ~ 0.3. The total luminosity 
change of RGB bump has been found to be ~1.83 in My (Ferraro et al. 1999) and ~1.76 
in Mk^ (Ferraro et al. 2000) in the same metalficity range. So, we derive the RGB bump 
luminosity variation in Mk^ with [Fe/H] as ^^f^. rll at^h] ~ ^-^^ metalficity 
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range [Fe/H] = —2.3 ~ 0.0, in the age range of 7 Gyr ~ 16 Gyr, and in the hehum abundance 
range y = 0.1 ~ 0.3. 

According to the theoretical luminosity function analysis of 46 RGB models of Sweigart 
& Gross (1978) by Yun & Lee (1979), when mass M = 0.9 Mq, the RGB bump bolometric 
luminosity variation with helium abundance Y is ^ ~ —5.0 in the helium abundance range 
Y — 0.1 ~ 0.3 and in the metallicity range [Fe/H] = —2.4 ~ —0.3. Here m is bolometric 
magnitude which is proportional to —2.5 x log(L/L0) and taken to be at the RGB tip. 
Comparison of the theoretical RGB bump bolometric variation rate by Yun & Lee (1979) 
with the theoretical RGB bump luminosity variation rate in My by Yi et al. (2001) gives 
Am f» 1.2lAMy ^ 1.26AMk,. So, when mass M = 0.9 Mq, the RGB bump luminosity 
variation in M^^ with helium abundance Y is found to be ^^y" ~ ~4.0 in the helium 
abundance range y = 0.1 ~ 0.3 and in the metaUicity range [Fe/H] = —2.4 ~ —0.3. 

According to Yi et al. (2001) and Ferraro et al. (1999), the RGB bump luminosity 
variation in My with age is fa 0.04 in the age range of 7 Gyr ~ 16 Gyr and in the 
metallicity range [Fe/H] = —2.3 ~ 0.0, where tg represents age in Gyr unit. Since the 
relation between the total luminosity changes of the RGB bump in My and in Mk^ is 
AMy 1.04AMx,, RGB bump luminosity variation in Mk^ with age is found to be 
^ 0.04 in the age range of 7 Gyr ~ 16 Gyr and in the metallicity range [Fe/H] = —2.3 ~ 
0.0. 

Applying equation (11) of Buzzoni et al. (1983) for the relation between R (A^hb/A^rgb, 
the ratio of the number of HB stars to the number of RGB stars brighter than the HB) and 
helium abundance Y to R values of 26 GGCs (Zoccali et al. 2000), the maximum range of 
helium abundance in GGCs at a given [Fe/H] is ~0.06 {AY = ±0.03). However, according 
to the Y data of Sandquist (2000) for 43 GGCs, the maximum range of helium abundance in 
GGCs at a given [Fe/H] is ~0.10 {AY = ±0.05). So, the maximum RGB bump luminosity 
variation in Mk^ caused by helium abundance variation in GGCs at a given [Fe/H] is found 
to be AMk, ~ -4.0Ay ^ ±0.20 mag. 

However, the maximum age spread of GGCs at a given [Fe/H] is ~6 Gyr {A{tQ) = ±3 
Gyr) according to Chaboyer et al. (1996), Richer et al. (1996), Salaris & Weiss (1997), 
Buonanno et al. (1998), Rosenberg et al. (1999), and VandenBerg (2000). So, the maximum 
RGB bump luminosity variation in Mk^ caused by age variation in GGCs at a given [Fe/H] 
is found to be AMk, ~ OMA{tg) ^ ±0.12 mag. 

Therefore, the maximum variation of RGB bump luminosity in M^^ with helium abun- 
dance and age spreads in GGCs at a given [Fe/H] is AMk, ~ ±\/0.202 ± 0.12^ ±0.24 
mag. This implies that the relations between the absolute magnitude of the RGB bump 
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(M^"™^) and metallicity [Fe/Hj^Qg^ or [M/H] derived in this work can be used to determine 
distance moduli of other GGCs with an uncertainty of approximately AM^:^ ^ ±0.24 mag 
at a given metallicity 



4. Summary 

Using photometric data from the 2MASS second incremental release point source cat- 
alog we found RGB bump features in 13 GGCs in Kg vs. (J — Kg) CMDs and measured 
accurate positions of RGB bumps for 11 GGCs excluding M69 and u Cen. Wc have found 
clear relations between M^"™^ and metallicity [Fe/H]cG97 or [M/H], thereby independently 
confirming the results of Ferraro et al. (2000). 

Combining the present sample and that of Ferraro et al. (2000), we extend the number 
of GGCs to 16 whose RGB bump positions have been measured, and determined more 
robust correlations between the absolute magnitudes of RGB bumps (M]^^'"^^) and metallicity 
[Fe/H]cG97 or [M/H]. Furthermore, these equations can be used to determine distance moduli 
of other GGCs with an uncertainty of AM^^ pa ±0.24 mag at a given metallicity. 

This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. This work is supported by the 
BK21 project of Korea through School of Earth and Environmental Sciences (Astronomy 
Program), Seoul National University. 
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Fig. 1. — Kg vs. (J — Ks) CMDs of GGCs whose RGB bump positions are accurately 
measured. 
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Fig. 1. — Continued. 
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Fig. 2. — Luminosity functions of GGCs in Fig. 1. In each diagram the upper part is the 
integral luminosity function and the lower part is the differential luminosity function. Arrow 
in each diagram indicates RGB bump position. The dashed lines in the upper part of each 
diagram are linear fits to the regions above and below the RGB bump. 
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Fig. 2. — Continued. 




Fig. 3.— (a) Ks vs. (J - Ks) CMD of M30. (6) K, vs. (J - Ks) CMD of M55. (c) The 
upper part is the integral luminosity function and the lower part is the differential luminosity 
function of M30. (d) Same as Fig. 3c but for M55. 
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Fig. 4. — (a) Relation between M^"™'' and [Fe/H]cG97 derived in this work. (6) Relation 
between M]^"™^ and [M /H] derived in this work. Solid line in each diagram is from eqs. (5a) 
and (5b). Small filled circle with error bar in the lower right corner of Fig. 46 shows typical 
error sizes of the values of small filled circles in Fig. 4. 
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Fig. 5. — (a) Relation between Af^^™'' and [Fe/H]cG97 combining results of this work and 
those of Ferraro et al. (2000) transformed into 2MASS system. (6) Relation between M^^™'^ 
and [M/H] combining results of this work and those of Ferraro et al. (2000) transformed 
into 2MASS system. Small open circles represent values from this work and small filled 
circles represent those of Ferraro et al. (2000) transformed into 2MASS system. Small 
filled triangles represent new values for NGC 6553 and NGC 6528 caused by mctallicity 
revision due to high resolution spectroscopic measurements by Carretta et al. (2001) and 
are connected to the original values of Ferraro et al. (2000) by dotted lines. Two small circles 
in one larger circle represent different values of same GGC common in both works and larger 
circle represents average of values of two small circles enclosed in the larger circle. In the 
case of Ml 07 the two values differ only by 0.01 mag so they are nearly indistinguishable. 
Solid line in each diagram is from eqs. (6a) and (6b) and dashed line in each diagram is from 
equations given in each panel of Fig. 13 of Ferraro et al. (2000) transformed into 2MASS 
system and dot-dashed hue in each diagram is from eqs. (5a) and (5b) as given in Fig. 4a 
and 4&. Small filled circle with error bar in the lower right corner of Fig. 5b shows typical 
error sizes of the values of small open and filled circles and small filled triangles in Fig. 5. 
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Table 1. Parameters of M22 and M2. 



M22 M2 



Parameter Value Reference Parameter Value Reference 



[Fe/H]ca97 


-1.48 


1 


[Fc/H]cG97 


-1.46 


5 


[M/H] 


-1.25 


2 


[M/H] 


-1.26 


2 


<^HB> 


14.10 


3 


<Vhb> 


15.93 


5 


VzAHB 


14.16 


2 


VzAHB 


15.99 


2 




0.59 


2 


^/^,^ZAHB 


0.59 


2 


E{B - V) 


0.34 


4 


E(5 - V) 


0.02 


5 


(m - M)CG97 


12.52 


2 


(m - M)CG97 


15.34 


2 



References. — (1) Carretta & Gratton (1997); (2) this work; (3) Peter- 
son & Cudworth (1994); (4) Harris (1996); (5) Lee & Carney (1999) 



Table 2. Parameters and derived RGB Bump Positions of Program GGCs. 



Cluster 


[fo/h: 




[M/ 


H] 


E(D 


- V) 


{rri — 












71 rBiirrip 




M15 (NGC 7078) 


-2, 


.12 


-1, 


.91 


0, 


.09 


15. 


,15 


12, 


.95 


± 


0, 


.05 


-2, 


.23 


± 


0, 


.21 


M22 (NGC 6656) 


-1, 


.48 


-1 


.25 


0, 


.34 


12. 


,52 


10, 


.80 


± 


0, 


.10 


-1 


.83 


± 


0, 


.22 


M2 (NGC 7089) 


-1, 


.46 


-1, 


.26 


0, 


.02 


15. 


,34 


13 


.55 


± 


0, 


.05 


-1 


.80 


± 


0, 


.21 


M13 (NGC 6205) 


-1 


.39 


-1 


.18 


0, 


.02 


14. 


,43 


12, 


.55 


± 


0, 


.05 


-1 


.89 


± 


0, 


.21 


M4 (NGC 6121) 


-1 


.19 


-0, 


.94 


0, 


.36 


11. 


,68 


10, 


.00 


± 


0, 


.10 


-1 


.80 


± 


0, 


.22 


NGC 362 


-1, 


.15 


-0, 


.99 


0, 


.05 


14. 


,68 


13 


.15 


± 


0, 


.05 


-1, 


.55 


± 


0, 


.21 


NGC 1851 


-1, 


.08 


-0, 


.88 


0, 


.02 


15, 


,46 


13, 


.80 


± 


0, 


.10 


-1, 


.67 


± 


0, 


.22 


NGC 288 


-1, 


.07 


-0, 


.85 


0, 


.03 


14. 


,73 


13, 


.25 


± 


0, 


.05 


-1, 


.49 


± 


0, 


.21 


M107 (NGC 6171) 


-0, 


.87 


-0, 


.70 


0, 


.33 


13, 


,95 


12, 


.50 


± 


0, 


.10 


-1, 


.56 


± 


0, 


,22 


47 Tuc (NGC 104) 


-0. 


,70 


-0, 


.59 


0, 


.04 


13, 


,32 


12, 


.10 


± 


0, 


.10 


-1, 


.23 


± 


0. 


,22 


M71 (NGC 6838) 


-0. 


.70 


-0, 


.49 


0, 


.25 


12, 


,97 


11 


.95 


± 


0, 


.05 


-1 


.10 


± 


0. 


,21 



Note. — Metallicities, interstellar reddenings, and distance moduli of GGCs are taken from Table 
1 of Ferraro et al. (1999) except for M22 and M2 whose values are taken from Table 1 of this work. 
More details are described in the text. 
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Table 3. Comparison of RGB Bump Positions of 3 overlapping GGCs in Ferraro et al. 

(2000) and in This Work. 





Ferraro et al 


. (2000) 


This work 


Average^ 


Cluster 


^Bump 


^Bump 


^Bump j^Bump 

S Kg 




M15 
M107 

47 Tuc 


12.86 ± 0.05 - 
12.51 ± 0.05 - 

12.01 ± 0.05 - 


2.32 ± 0.21 
1.55 ± 0.21 

1.32 ± 0.21 


12.95 ± 0.05 -2.23 ± 0.21 
12.50 ± 0.10 -1.56 ± 0.22 

12.10 ± 0.10 -1.23 ± 0.22 


-2.28 ± 0.15 
-1.55 ± 0.15 

-1.28 ± 0.15 



Note. — Ferraro et al.'s (2000) original results in Glass system were transformed into 2MASS 
system according to equation (Al) of Carpenter (2001). 

^Weighted mean average by errors. 
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Table 4. Absolute Magnitudes of RGB Bumps of 8 GGCs derived by Ferraro et al. 

(2000). 



Cluster 


[Fe/H]cG97 


[M/H] 


^Bump 


^Bump 


M15 (NGC 7078) 


-2.12 


-1.91 


-2.28 ± 0.21 


-2.32 ± 0.21 


M68 (NGC 4590) 


-1.99 


-1.81 


-2.30 ± 0.21 


-2.34 ± 0.21 


M55 (NGC 6809) 


-1.61 


-1.41 


-1.94 ± 0.21 


-1.98 ± 0.21 


M107 (NGC 6171) 


-0.87 


-0.70 


-1.51 ± 0.21 


-1.55 ± 0.21 


47 Tuc (NGC 104) 


-0.70 


-0.59 


-1.28 ± 0.21 


-1.32 ± 0.21 


M69 (NGC 6637) 


-0.68 


-0.55 


-1.15 ± 0.21 


-1.19 ± 0.21 


NGC 6553 


-0.06^ 


+0.08*^ 


-0.75 ± 0.22 


-0.79 ± 0.22 


NGC 6528 


+0.07^ 


+0.22^ 


-0.54 ± 0.22 


-0.58 ± 0.22 



Note. — Cols. (2) and (3) list metallicities taken from Table 1 of Ferraro 
et al. (2000) except for NGC 6553 and NGC 6528. Col. (4) lists absolute 
magnitudes of RGB bumps in Glass system and col. (5) lists those values 
transformed into 2MASS system according to eq. (Al) of Carpenter (2001). 

^These values were taken from Carretta et al. (2001). 

''These values were calculated by equation (3) in § 3.1 and more details 
are described in the text. 



